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Abstract 

By means of multiscale molecular simulations, we show that telechelic-star polymers are a simple, ro- 
bust and tunable system, which hierarchically self- assembles first into soft-patchy particles and then into 
targeted crystalline structures. The self-aggregating patchy behavior can be fully controlled by the num- 
ber of arms per star and by the fraction of attractive monomeric units at the free ends of the arms. Such 
self-assembled soft-patchy particles while forming, upon augmenting density, gel-like percolating networks 
and stable ordered structures, preserve properties as particle size, number and arrangement of patches per 
particle. In particular, we demonstrate that the flexibility inherent in the soft-patchy particles brings forward 
a novel mechanism that leads to the stabilisation of diamond and simple cubic crystals over a wide range of 
densities, and for molecular sizes ranging from about 10 nm up to the micrometer scale. 

PACS numbers: 82.70.-y, 81.16.Fg, 05.10.-a, 83.80.Uv 
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The current design of materials in the nanometer scale with its vast applications for, e.g., molec- 
ular electronics, therapeutic vectors and diagnostic imaging agent carriers, or photonics, focuses 
on the fabrication of a variety of organic and inorganic building blocks with different sizes and 
shapes, which are chemically decorated with discrete and specific interaction sites. Over the past 
years, functionalised particles gained an important role as building blocks for various materials 
Q3-0, rising the need for fabrication techniques that are both simple enough, and, at the same 
time, allow for a large-scale and reliable production of these particles. State of the art methods 
for the synthesis of such particles include lithography ASH, micro fluidics JE1 or glancing angle 
deposition flUEl. However, shape- specific particle fabrication becomes more challenging when 
the particle has to be decorated by a large number of patches or if the shape and location of the 
patches has to be tuned and controlled with extremely high precision. 

In this work, we show that di-block copolymer stars are a macromolecular system that allows 
for an extremely well-controlled production of novel particles with a potentially unlimited num- 
ber of patches. Diblock copolymer stars are star-polymers whose arms are characterized by dual 
physical or chemical functionalities. In particular, each arm is an amphiphlic A5-diblock with a 
solvophilic, athermal part A that is grafted in the center of the star, followed by a solvophobic, 
thermosensitive part B. Accordingly, we adopt the terminology telechelic star polymers (TSP) 
throughout. The solvophilic nature of the interior monomers results in a preferential maximiza- 
tion of the exposure to the solvent, whereas the solvophobic tail monomers tend to minimize their 
interaction with the same. By employing a coarse-grained approach, we surpass the molecular 
size limits of previous computational and theoretical approaches IfTTl [T2l demonstrating how, for 
extremely large and very dense systems (up to 10 7 monomers), by modifying the chemical com- 
position of the polymer chains and the number of arms, it is possible to obtain stable soft-patchy 
particles that survive in their self-assembled patchy state up to high concentrations. In contrast to 
the conventional hard-patchy colloids |[T3Tl . the patches here are formed by the self-association of 
flexible polymer chains, which can recombine and rearrange. The process of self-assembling in 
patches is then succeeded by a hierarchical assembly on an even larger scale. The properties of 
TSP's are determined by the interactions between the monomers and the solvent, the percentage a 
of solvophobic monomers per arm, and the number of arms /. In what follows, each of the / arms 
is composed by a total of Ma solvophilic and Mg solvophobic monomers. 

By using a first-principles coarse-graining procedure based on a soft-effective segment pic- 
ture [14] (see Supplemental Material II291D , we represent the TSP's by modeling the AB blocks of 



2 



every arm with ha and ng blobs. Each individual blob (or effective segment) contains a fraction 
niA = Ma/ua or m# = Mg/ng monomers respectively. As was shown in Refs. [fT4l IT3T1 . once the 
set of potentials characteristic for a specific AB diblock at a given temperature is computed, any 
asymmetry ratio a = Ma/ (Ma+Mb) can be obtained by adjusting the numbers ha and ng of blobs 
in each arm. Choosing the number of monomers in every blob large enough to be in the scaling 
regime enables us to obtain reliable properties and a general behavior of the polymeric particles, al- 
lowing thus for a mapping onto experimental systems. The underlying microscopic molecules an- 
alyzed here are formed by M = 10000 monomers per arm, which are represented in our approach 
by at least n m i n = 50 blobs. Additional simulations of molecules with M = 1000,2000,5000 and 
7000 monomers per arm have produced identical results, therefore implying the scalability of the 
properties of the system with the size of the molecule, as long as each arm has sufficient monomers 
to be in the scaling regime. To extract the set of effective potentials representing the effective seg- 
ments within the arms, full monomer simulations of diblock copolymers were performed at zero 
density employing a lattice model. In the full monomer representation, the solvophilic part was 
represented by self avoiding random walks, while to represent the solvophobic part we used a 
square-well interaction between monomers, adjusting the well depth so that the 5-blocks of the 
chains lie slightly below their ©-temperature. The potentials, extracted taking into account all the 
many-body interactions acting in the system 031[T6]|, 316 depicted in Fig.[T] 

We first characterize the system in the infinite dilution limit. To this end, the self- aggregation 
properties of a single molecule are classified as a function of its number of arms / and the per- 
centage a of attractive monomers per arm. Hereto, single molecule simulations are performed 
within the coarse-grained approach described above for various combinations of / G [2, 20] and 
a G [0.1,0.8]. The single-molecule conformational diagram showing the soft patchy particle na- 
ture of the TSP's and the tunability of the number of patches p as a function of / and a is presented 
in Fig. [2j It reveals that molecules with a < 0.3 show no sign of self-aggregation, whereas such 
a behavior can clearly be identified for a > 0.3. For given (/, a) values, molecules self-assemble 
into patchy particles characterized by as many as six patches in the parameter range considered. 
The number of patches that are formed and the concomitant anisotropy of the particles, can be 
precisely tuned by choosing an appropriate combination of (/, oc). Previous full monomer works 
on TSP's with much shorter arms (M ^ 30) IITTl[T7ll showed that, for a certain range of the (/,0c) 
parameters, these molecules exhibit similar behavior as the stars with M = 10000 monomers per 
arm: they self-assemble into structures with single or multiple assembly sites that were named 
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Figure 1: Color online. The non-bonded potentials Vaa in grey, Vab in blue and Vbb in green and the 
tethering potentials cp^A in grey, ($>ab in blue and ($>bb in green between the various blobs of type A and 
B. Color coding of the blobs: gray (solvophilic), green (solvophobic), anchor point (red). The radius of 
gyration of the A and B blobs is the same and labeled as r g . 

single or multiple "watermelons". The full monomer analysis showed that the single collapsed 
watermelon conformation becomes unstable for short chains at fixed temperature when increasing 
the number of arms (/ > 6) in favor of multiple watermelon structures, identical to the patchy 
particles found here. 

Recently, it has been demonstrated that linear block copolymers have the ability to assemble 
into micelles or cylinders decorated by patches lfT8l[T9l . In Ref. [fT9ll . it was shown that an increase 
of the density of block copolymers in solution induces a growth in the patchy aggregate, whose 
size and patchiness therefore change with concentration. If one aims at producing robust soft 
patchy particles, however, the stability of the aggregate formed by block copolymers should be 
tested against the density of copolymers in solution. Therefore, to assess to which extent TSP's are 
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Figure 2: Single-molecule conformational state diagram of TSP's, shown as a function of the number of 
arms / and asymmetry ratio a. The number of multi-arm patches p is indicated on the plot. Stars with 
a < 0.3 remain in an open configuration. 

stable, self-aggregated patchy particles that maintain their zero density properties even in the semi- 
dilute and concentrated regimes, simulations at finite density have been performed for a number of 
selected (/, a) -combinations. We spaced our simulations over a wide range of densities, starting 
from 0.5p* to 10p*, where p* = 3/(4%Rg) is the overlap density of TSP's having gyration radius 
R g . Simulations ofNe [100, 250] stars made of / G [3, 20] arms that have a percentage of attractive 
monomers a 6 [0.2,0.8] have been performed. 

For every density, we quantify the properties of the single particle and of the whole solution by 
performing both a cluster analysis and a topological analysis by means of the Euler characteristic. 
Both approaches (see the Supplemental Material for details on the methodologies) consistently 
reveal that the finite density aggregation phenomenon is a consequence of the structure of the 
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Figure 3: Snapshot of a percolating gel formed by telechelic star polymers with / = 5 arms and 60% of 
attractive monomers per arm. Color coding as in Fig. [T] The snapshot shows the simulation box and some 
parts of the periodically repeated system. 

local intra-particle aggregation that leads to the formation of patches. As concentration grows, 
inter-star patch merging leads to the formation of low density percolating gels or networks, see 
Fig. [3] The cluster analysis of the gel shows that the zero density single star properties, i.e., the 
average number of patches per star, the average angle between the patches, and the distance of 
the patches from the centre of the star, together with the fluctuations over the average values, are 
unaffected by the density of stars in solution. This is further supported by the topological analysis. 
Here, covering spheres of varying radius R are placed at the center of each solvophobic blob 
and the Euler characteristic of the resulting pattern, %(R), is considered as a function of R. Typical 
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results are shown in Fig.S2 in the Supplemental Material. The salient feature that can be seen there 
is the existence of broad, flat plateaus in the range 10 < R/r g ^ 50, which are the signature of the 
existence of well-formed and well-separated patches. The number of patches per TSP, the number 
of blobs per patch as well as the average distance between the patches can all be extracted from 
this topological invariant and the results show a fully consistent picture with that obtained at the 
infinite dilution limit. Thus, even in the random gel, the zero-density self-aggregation properties 
of TSP's remain unaffected. 

TSP's with a given patchiness thus emerge as novel candidates for the stabilization of ordered, 
crystalline states possessing local coordination compatible with the given particle anisotropy. We 
therefore consider two target crystals to assess the large scale self-assembling properties of TSP 
soft patchy particles, namely a diamond and a simple cubic crystal. Fourfold-coordinated TSP's 
are compatible with a diamond lattice, whereas sixfold-coordinated stars are a natural choice to 
obtain stable simple cubic crystals. Hard core patchy particles have proven to be an extremely 
versatile system for the synthesis of a wide variety of crystals; however even if several geometries 
of patches have been considered, some crystals have been found to be stable only within a very 
small interval of the parameters range, i.e., number of patches, strength of the attractive potential 
and temperature [|4ll20l. 

We found that the flexibility both in the stiffness of the core and in the localization of the 
patches on the particle is a major advantage when patches undergo large-scale self-assembly 
processes. This flexibility has its main origin at the entropically driven fluctuations of the self- 
avoiding blocks of the arms. In particular, we placed the anchoring points of the TSP's at the 
positions of various perfect lattices (simple cubic, diamond, fee, bec) and their arms at various 
different initial equilibrated configurations and we performed very long simulations, letting the 
system free to equilibrate and evolve. Since the diamond lattice has a coordination of four neigh- 
bors, the natural choice for the (/, a) parameters of the constituent stars, are / in the range [10, 16] 
and a £ [0.4,0.8], values for which, at zero density, stars self assemble in four patches. Similarly, 
the simple cubic requires stars with / = 20 and a £ [0.5,0.7], that self assemble in six patches 
at low density. However, in principle, stars could stabilize different crystal structures whose co- 
ordination number is compatible with the number of arms per star instead of with the number 
of patches that are formed at zero density. If the stabilization mechanism is dominated by the 
attractive monomers at the end of every arm per star, for example for the range of / that assem- 
bles in four patches, then indeed the diamond, bec and fee crystals (with coordination numbers 
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of 4, 8 and 12 respectively), could be all viable candidates. To reduce the computational cost, we 
considered the scenarios corresponding to / = 15 and a = 0.6 and / = 15 and a = 0.4 for the 
fourfold-coordinated TSP's, and / = 20, a = 0.5 for the sixfold-coordinated ones at a wide range 
of reduced densities p/p* G [1.3, 10] . 

For all the densities analysed, both bcc and fee crystals melted within the simulation time; on 
the contrary, the diamond and simple cubic crystals remain mechanically stable respectively for 
the four- and sixfold-coordinated particles. A mean square displacement (MSD) analysis of the 
crystals (not shown) has been performed to assess the stability of the system; both the bcc and fee 
displayed a continuously growing MSD, while the diamond and simple cubic crystals analysed 
show saturation to a plateau at long times. A cluster analysis performed on the crystal shows that 
the zero density single-particle properties are preserved in the crystalline phase for a wide range 
of densities. In particular, the average number of patches (patchiness, p), the bond angles between 
the segments connecting the anchoring point to two different patches (inter-patch angle (0) as well 
as the distance from the anchoring point to the center of mass of the patch (patch extension L) were 
fully preserved in the crystalline phase for a wide range of densities with respect to the values they 
have at infinite dilution, see Tables ST1 and ST2 of the Supplemental Material. 

Tetrahedrally coordinated TSP's, such as the ones with / = 15 considered here, stabilize di- 
amond crystals by merging the attractive patches belonging to molecules on neighboring sites, 
and sixfold-coordinated TSP's, such as the ones with / = 20 considered, stabilize simple cubical 
crystals via a similar mechanism. The flexibility in the location of the patches and in the distance 
between the patches and the anchor point induces a novel mechanism for the stabilization of the 
two crystals. Patches are not only merging among the first neighbors but involve further stars all 
the way to the third neighbors. Such a preference can be understood by taking into account the 
fact that the TSP's maintain their patchy structures at all crystal densities, hence with decreasing 
lattice spacing the stars are forced to pack by pushing the bulky self-avoiding part of the arms 
into the empty spaces around each lattice site. For the diamond crystal, such spaces are distributed 
with tetrahedral symmetry, compatible with the natural angular arrangement of the patches formed 
by each star. Once pushed into such spaces, the attractive patches can interact with the patches 
belonging to the other stars present in the same volume, which correspond up to the third neigh- 
bors. Similarly to the diamond, the simple cubic crystal is also stabilised by the packing of the 
self-avoiding branches of the TSP's, which reach patches in the neighboring cells. Snapshots of a 
typical diamond crystal formed for (/,0c) = (15,0.4) and of a stable simple cubic crystal formed 
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Figure 4: Stable diamond lattice (a) and simple cubic crystal (b) obtained respectively for the / = 15 , a = 0.4 
and / = 20, a = 0.5 systems. The red lines connecting the anchor points of the stars, are a guide to the eye, 
drawn to stress the four fold coordinated diamond lattice structure and six fold coordinated simple cubic 
structure. Color coding as in Fig.[T](Bonds not shown). 

by (/, oc) = (20,0.5) are shown in Fig. |4j In the figure, the diamond and the simple cubic ar- 
rangements of the centres of the stars, surrounded by a high density cloud of polymers, are clearly 
visible. 

The same mechanism is also responsible for the rapid melting of the bcc and fee crystals which 
are not compatible neither with the packing of the self-avoiding arms and lack the ability to create 
local environments that favor a tetrahedral arrangement. This effect can be observed by monitoring 
the preference for first second and third neighbor patch fusion at increasing densities. Moreover, 
and in contrast to the case of athermal star polymers and chains, the radius of gyration of the 
star is only very weakly affected by increasing the concentration. In fact, even at the infinite 
dilution the size of the multi-patch TSP's considered here is significantly smaller than that or their 
athermal counterparts with the same / and number of monomers per arm. Concomitantly, no 
further shrinkage is necessary as the density grows and the TSP's maintain the size and form they 
have at the infinite dilution limit. 

We have shown that telechelic star polymers are extremely tunable and versatile building blocks 
that have the ability to hierarchically self-assemble into soft patchy particles, which organize them- 
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selves into complex structures. Using advanced computational methods we fully characterize the 
self-organization properties of telechelic star polymers from zero density all the way up to high 
concentration regimes, where we observe stable diamond and simple cubic crystals over a wide 
range of densities. At low density the stars self-aggregate into particles decorated by a number of 
patches that is completely tunable by choosing different combinations of the functionality / of the 
stars and on the percentage a of attractive monomers per arm. Results remain invariant within the 
hmits M E [10 3 , 10 4 ] monomers per star, implying their their scalability with molecular size. In 
addition, the major findings of our work rest on the competition between the enthalpic attractions 
between the solvophobic, terminal blocks and the self-avoiding central ones. Therefore, we expect 
that our results will remain essentially unaltered over a substantial range of temperatures T < T®, 
provided the terminal monomers do not reach the fully collapsed regime. 

The desired patchiness naturally arises as the result of self- organization of the particles at the 
molecular level, rendering laborious preparation techniques unnecessary. These self-assembled 
patchy structures survive at intermediate density regimes, where the particles form percolating 
networks. Telechelic stars therefore have enormous potential as simple, extremely tunable and 
controllable soft and flexible patchy particles. Given their stability at finite density, such parti- 
cles are excellent candidates as tunable building block for new materials and they provide both 
selectivity in the crystal structures they stabilize and scalability in the lattice constant, since they 
maintain their properties within the range of a few dozens to hundreds of thousands of monomers 
per arm. In contrast to hard patchy colloids, soft patchy particles take advantage of entropy in 
selecting and stabilizing structures. In particular, allowing for fluctuations in the localization of 
the patches, results in an increase in the entropy of the target structures therefore enhancing their 
stability over a wide range of densities. 

State of the art of technology allows to fully control the synthesis of telechelic star polymers, 
and various experiments have been performed on solutions of telechelic star polymers at finite 
density ET1I221 showing an extremely interesting aggregation scenario. With this study we have 
established telechelic star polymers as a novel, promising, and extremely versatile self-assembled 
and tunable patchy particle system. The capability of these macromolecules to self-assemble into 
robust, soft patchy particles that maintain their shape and directionality deep into the semidilute 
regime, surpasses problems related to the fabrication of conventional patchy colloids, opening-up 
ample new opportunities for experiments and technological applications. Moreover, their flexibil- 
ity, which stems from their soft nature as fluctuating polymers, selectively stabilizes target crystals 
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for a broad range of concentrations, opening the way for the self-assembly of a variety of desired 
solids. 
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SUPPLEMENTAL MATERIAL FOR PAPER: "TELECHELIC STAR POLYMERS AS SELF- 
ASSEMBLING UNITS FROM THE MOLECULAR TO THE MACROSCOPIC SCALE" 



Coarse-graining method 

The coarse graining used throughout this work, is an extension of the soft effective segment 
representation (SES) previously developed for both homopolymer linear (grafted and ungrafted) 
chains Il23l 1241 and block copolymers Ifl4l . The methodology has been proven extremely effi- 
cient in reproducing fundamental physical properties of polymers in the semi-dilute solution. The 
strategy consists in representing the polymer as a chain of effective segments or blobs. Each 
blob, representing a group of monomers, interacts with other blobs via effective interactions that 
have been computed for dimers (at infinite dilution) taking into account many-body contributions 
lfl4l[T6"1l . The effective interaction computed at zero density successfully reproduces the physical 
properties of the system (e.g. pair distribution functions) till the blobs start to overlap, e.g., till the 
density p/, of blobs in solution reaches the overlap density of blobs p£. The SES strategy consists 
in representing the polymer chains with a number of blobs high enough to ensure that the condition 
pb < Pi is satisfied. Hence, the higher the density in solution the finer is the level of the coarse 
graining required. In the range of densities considered in this paper we found that the minimum 
number n m { n of blobs that is required for obtaining a proper description in the scaling limit of 
telechelic star polymers, is determined by the local monomer density within the star. Therefore, to 
tune the level of coarse graining needed to reproduce the physical quantities in the system while 
guaranteeing that the coarse graining is not affecting the results, we performed an extensive anal- 
ysis of single star properties, upon increasing the number of blobs per arm, therefore augmenting 
the detail of the coarse graining that we are using. 

We found that for n > n m { n all the properties of the single star (radius of gyration of the star, 
radius of gyration of the arms, radius of gyration of the two components of the arms, intra- 
aggregation properties), are independent of the level of the coarse graining used to represent the 
star, thus defining n m [ n as the minimum number of blobs needed for the results to be independent 



on the level of coarse graining used. We show in Fig. SI the estimate of the various radii of gyra- 
tion calculated with four different level of coarse graining, namely 50, 55, 60 and 65 blobs per arm 
for stars with / G [3, 9] and a = 0.4. As evident from the figure, the four different coarse graining 
levels perfectly agree in reproducing the radius of gyration of the whole star, the radius of gyration 
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Figure SI: Radius of gyration of whole star (black triangles) of each arm (red diamonds) of the self 
avoiding head of each arm (blue squares) and of the attractive tail of the arms (green circles) obtained using 
50, 55, 60 and 65 blobs per arm for a star with / £ [3,9] and a = 0.4. The various identical symbols 
correspond to different number of blobs. The inset shows the comparison between the radius of gyration of 
the self-avoiding part of the polymer arms (red a = 0.8, black a = 0.6, blue a = 0.4) and results obtained 
with full monomer simulations for self-avoiding homopolymers of equivalent length (dashed lines). The full 
monomers results have been rescaled in terms of the SES representation using the scaling law introduced in 

of each arm of the star, of the self-avoiding head and of the attractive tail of each arm. Therefore, 
as long as the number of blobs n is bigger than the minimum number of blobs n m { n (found to be 
^min = 50 for the cases at hand), all the properties of stars and system become coarse graining 
independent. Moreover, it has already been shown for the case of grafted homopolymer brushes 
that the method of coarse graining allows for a controlled back-mapping onto full monomer results 
Il24l . In this work, we always used a total number of segments or blobs n > n m { n , thus satisfying 
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the constraint for the validity of this approach. Performing simulations for systems slightly below 
the 0-temperature for the solvophobic tails, and in the athermal solvent regime for the solvophilic 
heads, grants an attractive, effective potential between the B and a repulsive, Gaussian potential 
between the A blobs. The microscopical asymmetry ratio can also be rewritten in terms of the 
coarse grained model, which reads: 

n A (r gA ) l/VA 

(X — — — (1) 

n A (r gA ) l/VA +n B (r gB )^ 

where r gA and r g s are the radii of gyration, and v A and Vb are the Flory exponents of the A and B 
effective segments respectively. Here we chose r gA = r g g = r g . Note that equation Q assumes that 
the number of monomers in every segment is large enough to be in the scaling regime of polymers 
(e.g., large number of monomers per arm). 

The Euler characteristic 

Systems assembling in low density percolating networks (or gels) are characterized by an inho- 
mogeneous distribution of particles in space; such inhomogeneity makes gels an extremely com- 
plex system to characterize, since high order correlations play an important role. For this reason, 
a structural analysis of network forming systems are not straightforward: common characteriza- 
tions, essentially based on two particle correlations, do not give a full insight on the assembled 
status of the system and more complex analysis are required for such phases. 

A robust indicator that has been used to describe morphologies on vast scales, ranging from the 
large-scale structure of the universe [25 1 to the microscopic scale of molecules and atoms [|26ll27ll 
is based on the Euler-characteristic %. The latter is a topological invariant that was introduced by 
Leonhard Euler to characterize surfaces of polyhedra by counting their number of corners (ver- 
tices), edges and faces. For the topological characterization of the gels, we consider the surface 
that is formed by a collection of spheres with radius R centered on the position of each solvopho- 
bic blob in the TSP. The Euler characteristic %(R) that we obtain within our analysis is therefore 
an indicator for a family of surfaces formally parametrized by R. The complexity of the anal- 
ysed surfaces entails an enhanced sensitivity of this measurement with respect to many-particle 
correlations and provides insights in both local and global structure of the gel phase. Examples 
of appropriately normalized %(i?) at finite densities, only considering the solvophobic blobs, are 



given in Fig. [S2]of the main text for systems of telechelic star polymers with / = 3,5,7 and 10 
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Figure S2: The Euler characteristic %{R) for a system of telechelic star polymers with / = 3,5,7 and 10 
arms and asymmetries a = 0.4,0.6, and 0.8 at finite densities. The x axis is scaled with the radius r g of 
gyration of the blob. Panel (b) is a zoom of panel (a) in the region of the plateaus of the Euler characteristic. 

arms and a = 0.4,0.6 and 0.8. 

A change in the value of %(R) is brought about by a topological difference in the surface induced 
by a small increase in the radius R of the covering spheres. The simplest of the events causing 
those changes is found when two unconnected spheres start touching on increasing the value of 
R. In this particular case the topology is modified, because two, initially disjoint, surfaces have 
merged to a single one. In the same fashion, a plateau in %(R) indicates the absence of such 
events, which implies that at those length-scales "nothing" is happening, or that there is a balance 
between different types of events in the more clangorous range of distance scales. For a more 
complete description of the methodology we refer the reader to Refs. E71l28ll . 



The plateaus of %(R) in Fig. S2 in the range 10 < R/ r g ~ 50 are indicative of the formation 



of well-separated patches by the telechelic stars. The beginning of each plateau is related to the 
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distances between solvophobic blobs within the same patch, whereas the end and the height of a 
plateau correspond to the minimal distance between distinct patches and the average number of 
solvophobic blobs in a patch respectively. 

Cluster analysis 

The hierarchical self-assembling properties of TSP rely on the ability of the stars to keep their 
low-density assembled patchy particle structure upon augmenting densities. In the infinite dilution 
regime, the identification of the patches formed by a single star, together with the characterization 
of their angular and spatial distribution around the centre of the star, is straightforward. At finite 
densities on the other hand, particles start aggregating with one another and the more they interact, 
the more difficult it becomes to distinguish the contribution of each star to the aggregates as well 
as to characterise the single star self-aggregating properties. The Euler characteristic is capable 
of identifying the presence of clusters and their sizes, however it does not allow for a proper 
characterization of the patchy nature of each star. 

A simple visual analysis of simulations snapshots shows that the clusters formed by the solvo- 
phobic arm ends are nicely defined and isolated, suggesting that a cluster algorithm should be 
capable of precisely identifying the clusters. Fundamental properties are then measured in details, 
such as the cluster (patch) size and arrangement around each stars, as well as the counting and 
classification of the clusters resulting from the merging of distinct patches. The cluster analysis is 
then performed both at zero density and at finite density to fully characterize the properties both 
of the single star and of the solution. Clusters are defined by introducing a cut-off radius R cut : 
particles are classified as neighbours if their distance is below R cut . 

Clusters are then built by grouping particles that share at least one neighbour. According to 
this definition, a patch is a cluster of solvophobic particles. i? cut is a free parameter, that is chosen 
in a range of values around R cut ~ r g , where r g is the radius of gyration of the blob. The optimal 
i? cu t is such that the number and size of the clusters is not affected by a small variation AR cut of 
the cut-off radius. Such a cut-off value exists because the clusters are well separated in solution 
as it can be seen in the simulations snapshots (see Fig. 3 in the main article) and by the presence 
of plateau region in the Euler characteristic (see Fig. |S2| ). Once a list of cluster is obtained for 
each simulation snapshot, we measure the number of patches that belong to each star, their geo- 
metrical arrangement around the center of the star, and the degree of connectivity present in the 
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system. Such quantities are then averaged over the configurations sampled for each density and 
star parameters (functionality / and arm solvophobicity a). Accordingly, we can classify how the 
patches distribute around the stars and how they are shared among the stars in solution. The aver- 



ages are then compared to the zero density simulations as showed in Table ST1 and in Table ST2 



below, verifying the robustness of the patchiness p across a wide range of concentrations both in 
the ordered crystalline and in the disordered gel phase. 
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Table ST1: Properties of the TSP's in the diamond lattice for two different tetrahedrally coordinated 
molecules, namely / = 15, a = 0.6 and / = 15, a = 0.4 and for a six fold coordinated molecules, namely 
/ = 20, a = 0.5, over a wide range of densities, indicated at the first column. In the tables, p represents the 
average number of patches, co the bond angles between the segments connecting the anchoring point to two 
different patches, L the patch extension (the distance from the anchoring point to the center of mass of the 
patch) and at the last column, a denotes the bond length in the underlying monomer resolved system. 



/= I5,<x = 0.6 


/= 15,a = 0.4 


P/P* 


P 


CO [in degrees] 


L [in units of a] 


P/P* 


P 


CO [in degrees] 


L [in units of a] 


0.0 


4.01 


103.5 ± 27.8 


142.2 ± 29.2 


0.0 


4.02 


104.4 ± 31.8 


161.6 ±27.6 


1.3 


3.85 


102.6± 24.8 


135.4± 37.6 


1.9 


4.20 


100.7± 37.9 


165.0 ± 42.0 


2.2 


3.84 


102.6± 34.8 


132.2± 35.5 


3.0 


4.25 


101. 8± 38.0 


164.9 ±41.0 


3.8 


3.81 


103.5± 33.5 


132.8± 34.5 


5.2 


4.12 


102.7± 35.2 


168.9 ±38 


7.4 


3.86 


103.5± 34.3 


134.8± 32.2 


10 


4.13 


101. 8± 36.3 


168.1± 34.2 







/ = 20,a = 


0.5 


P/P* 


P 


co [in degrees] 


L [in units of a] 


0.0 


6.1 


96.5 ± 34.2 


168.0 ± 32.3 


1.28 


5.7 


94.1 ±38.4 


171.5 ±36.8 


3.05 


5.6 


93.4 ± 37.4 


172.7 ± 35.7 


5.72 


5.5 


94.6 ± 38.7 


169.4 ± 33.1 


10.3 


5.6 


88.4 ± 42.7 


155.5 ±27.1 
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Table ST2: Intra particle aggregation properties for star polymers with / € [5, 15] and fixed asymmetry 
a = 0.4, over a wide range of densities (infinite dilution, p/p* = 0, and in the gel phase, p/p* > 0). 



/= 


= 5 


/ = 


= 7 


/ = 


10 


/ = 


15 


p/p* 


P 


P/P* 


P 


P/P* 


P 


P/P* 


P 





1.12 





2.13 





2.84 





4.02 


0.10 


1.24 


0.25 


1.93 


0.38 


3.16 


0.52 


3.90 


0.20 


1.20 


0.50 


1.82 


0.75 


2.65 


1.03 


3.84 


0.25 


1.13 


0.63 


1.80 


0.93 


2.61 


1.30 


3.88 


0.30 


1.30 


0.76 


1.86 


1.12 


2.41 


1.55 


3.86 
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